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ABSTRACT 

Many  military  systems  must  be  capable  of  sustained 
operation  in  the  face  of  mechanical  shocks  due  to  projectile  or 
other  impacts.  The  most  widely  used  method  of  quantifying  a 
system’s  vibratory  transient  response  to  shock  loading  is  called 
the  shock  response  spectrum  (SRS).  The  system  response  for 
which  the  SRS  is  to  be  determined  can  be  due,  physically, 
either  to  a  collocated  or  to  a  noncollocated  shock  loading. 
Taking  into  account  both  possibilities,  one  can  define  the  SRS 
as  follows:  the  SRS  presents  graphically  the  maximum  transient 
response  (output)  of  an  imaginary  ideal  mass-spring-damper 
system  at  one  point  on  a  flexible  structure,  to  a  particular 
mechanical  shock  (input)  applied  to  an  arbitrary  (perhaps 
noncollocated)  point  on  the  structure,  as  a  function  of  the 
natural  frequency  of  the  imaginary  mass-spring-damper  system. 
For  a  response  point  sufficiently  distant  from  the  impact  area, 
many  Army  platforms  (such  as  vehicles)  can  be  accurately 
treated  as  linear  systems  with  proportional  damping.  In  such 
cases  the  output  due  to  an  impulsive  mechanical-shock  input 
can  be  decomposed  into  exponentially  decaying  sinusoidal 
components,  using  normal-mode  orthogonalization.  Given  a 
shock-induced  loading  comprising  such  components,  this  paper 
provides  analytical  expressions  for  the  various  common  SRS 
forms.  The  analytical  approach  to  SRS-determination  can  serve 
as  a  verification  of,  or  an  alternative  to,  the  numerical 
approaches  in  current  use  for  such  systems.  No  numerical 
convolution  is  required,  because  the  convolution  integrals  have 
already  been  accomplished  analytically  (and  exactly),  with  the 
results  incorporated  into  the  algebraic  expressions  for  the 
respective  SRS  forms. 

INTRODUCTION 

Modern  warfare  calls  for  many  military  systems  to  be 
capable  of  sustained  operation  under  extreme  environmental 
conditions.  The  mechanical  shocks  from  such  sources  as  blast- 
waves  and  projectile  impacts,  and  even  vehicular  motion  over 
rough  terrain,  make  high  demands  on  military  equipment.  For 
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design  and  analysis  purposes,  the  vibratory,  transient  response 
of  systems  (or  of  system  models)  subject  to  mechanical  shock 
is  typically  captured  using  two  frequency-dependent 
(“spectral”)  tools  [1,  2]:  (1)  the  Fourier  spectrum;  and  (2)  the 
shock  response  spectrum  (SRS). 

A  two-dimensional  SRS  [1]  (typically  termed  simply  an 
SRS)  represents  graphically  the  frequency  content  of  a 
specified  shock  input  d{t)  in  terms  of  the  maximum  response 
x(t)  it  would  induce  in  a  hypothetical,  single-degree-of- 
freedom  (SDOF)  mass-spring-damper  (MSD)  system, 
seismically  subjected  to  the  shock.  (Refer  to  Figure  1.)  The 
SRS  plots  a  selected  kinematic  measure  of  the  maximum  time- 
domain  motion-response  (of  mass  m )  against  the  SDOF- 
system  natural  frequency,  with  the  frequency  varied  over  some 


Figure  1.  Hypothetical  SDOF  MSD  system, 
for  SRS  determination 

range  of  interest.  In  producing  the  plot,  both  the  shock  input 
and  the  SDOF-system  damping  ratio  are  held  constant  across 
the  entire  range  of  natural  frequencies.  The  shock  input  is 
typically  a  position  input  [r/(f)  in  Figure  l] .  Some  useful 
measures  of  the  response  include  (1)  relative  displacement, 
x(t)-d{t)\  (2)  “pseudovelocity,”  a>„  [x(r)-r/(f)];  and  (3) 
absolute  acceleration,  x(t) .  As  used  in  this  report,  the 
respective  SRS’s  are  designated  and  defined  as  follows  [2]: 
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Spectral  displacement  (or  relative-displacement  SRS): 

■S’ d («„):=  •  (!) 

Spectral  velocity  (or  pseudovelocity  SRS)  [3,  4]: 

*SV  H )  :=  (»„  -  40  max = )  •  (2) 

Spectral  acceleration  (or  absolute  acceleration  SRS): 

SaM’-=  |40Lx-  (3) 

Consider  now  the  case  of  a  mechanical  shock  c(t )  applied 
to  an  arbitrary  point  C  of  a  generic  system  S.  (Refer  to  Figure 
2.)  If  a  nonlinear  finite-element  model  exists  for  S,  the  induced 
displacement  d(t)  at  some  other  point  D  on  S  can  typically  be 
determined  numerically,  and  then  used  to  calculate  (again, 
numerically)  the  indicated  spectral  quantities  at  D.  In  general, 
though,  analytical  evaluation  of  SD  ,  Sy  ,  or  is  not  possible 
for  nonlinear  systems.  However,  for  a  response  point 
sufficiently  distant  from  the  impact  area,  many  army  platforms 
(e.g.,  vehicles)  can  be  accurately  treated  as  linear  systems  with 
proportional  damping.  In  such  cases  analytical  SRS 
determination  proves  possible. 

Physical  Point  of  SRS 


Figure  2.  Shock-loaded  system  S, 
with  (hypothetical)  attached  MSD  system 
for  SRS  determination 


Consider  the  typical  case  of  an  undamped  or  underdamped, 
linear  system,  subjected  to  a  mechanical-shock  disturbance,  at 
point  C.  In  such  cases,  analysis  by  the  method  of  modal 
superposition  yields  an  expression  for  the  induced  displacement 
d(t)  (at  D )  as  a  linear  combination  of  modal  coordinates  [5,  6]. 
These  coordinates  are  the  solutions  to  their  respective 
decoupled  scalar  differential  equations.  The  inputs  to  these 
scalar  equations  are  modal  forces,  which  are  themselves  scalar 
multiples  of  the  mechanical-shock  disturbance.  If  the 
mechanical  shock  can  be  approximated  as  an  impulse,  the 
induced  displacement  d{i)  will  be  a  linear  combination  of 
damped  sinusoids.  Denoting  the  fth  sinusoid  by  d ,  ( t ) ,  one  can 
express  the  displacement  at  D  in  the  following  form: 

<*(')= (4) 

i= 1 

where  d,  (t )  =  Dte  a,‘  sin(o;f  +  <t>i)ujt).  (5) 


This  displacement  d(t)  serves  as  the  input  to  the  SDOF 
MSD  system  (see  Figure  2).  For  the  ;th  component  c/,  (/  ) .  the 
constants D,,  a,  ,  co, ,  and  (]),  represent  the  input  amplitude, 
rate  of  exponential  decay,  oscillation  frequency,  and  phase 
shift,  respectively.  Each  component  begins  at  time  t  =  0  (the 
time  of  impact),  as  indicated  by  the  unit  step  function  uA(t) . 
As  documented  below,  if  one  assumes  the  induced 
displacement  d(t)  to  be  given  by  Equations  (4)  and  (5),  one 
can  find  an  analytical  expression  for  the  relative  displacement 

5(f)  =  x(f)-<i(f),  (6) 

and  for  the  absolute  acceleration  x(t) ,  to  use  in  evaluating 
SD(a>n),  Sr (ff>„),  and  ^(co,,)  for  any  q  in  the  range 
0<q<l.  These  analytical  expressions  can  either  be  used  to 
plot  the  respective  spectral  quantities,  or  to  check  the  plots 
found  using  alternative  evaluation  methods. 

PROBLEM  STATEMENT 

Consider  the  linear,  SDOF  MSD  system  shown  above  in 
Figure  2.  Assume  the  displacement  d(t)  at  point  D  to 
comprise  a  linear  combination  of  V  exponentially  decaying 
sinusoidal  inputs  (i.e.,  to  the  SDOF  MSD  system),  as  given  by 
Equations  (4)  and  (5).  The  objectives  of  this  research  effort 
are  to  develop  analytical  forms  for  the  spectral  quantities 
S^(ff>„),  and  S^ro,,),  as  defined  by  Equations  (1), 
(2),  and  (3),  respectively.  The  general  case  (  0  <  q  <  1 )  will  be 
considered  first,  and  then  the  special  case  of  q  =  0  . 

SOLUTION  FOR  THE  GENERAL  CASE  ( 0  <  q  <  1 ) 

Basic  approach 

The  system  differential  equation  of  motion  (DEOM)  can  be 
expressed  by 

mx  +  cx  +  kx  =  cd  +  kd  .  (7) 

In  standard  form  (i.e.,  in  terms  of  the  damping  ratio  c;  and  the 
natural  frequency  CO/( ),  the  DEOM  is 

x  +  2q(Qnx  +  (o2nx  =  2q(ond  +  (o2nd .  (8) 

Since  d(t )  is  known  [Eqs.  (4),  (5)],  one  can  find  the  relative 
displacement  by  solving  Equation  (8)  for  x(t) .  Then  Equations 
(1),  (2),  and  (3)  can  be  used  to  formulate  the  respective  shock 
response  spectra. 

Response  x(t)  to  a  generic  induced  displacement  d(t ) 

Consider  first  the  response  x(t)  of  the  SDOF  MSD  system 
to  an  arbitrary  (assumed  Laplace-transformable)  induced 
disturbance  d(t ).  Let  the  Laplace  transforms  of  x(t)  and  d(t) 
be  represented  by  X{s)  and  D[s) ,  respectively.  Then  the 
Laplace  transform  of  Equation  (8)  is 

s2X(s )  -  +) -  +)  +  2  qro„  [Of(.sj  -  x(0  +)]  +  co  2,X(s)  ^ 

=  2  qcon  [.s’ £>(?)  -  d(0  +)]  +  . 
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Rearrangement  yields 


40- 


2q(ans  +  cd- 
s2  +  2qans  +  ft)2 


40- 


s  +  2qft>  „ 


s2  +  2qans  +  k>2 


s2  +  2qe>Hs  +  o)2 


i(()  +)  - 


2?ro„ 


s2  +  2qo jns  +  ft)2 


r(0+) 

d  (O  + ) . 


(10) 


Taking  the  inverse  Laplace  transform  of  Equation  (10),  one 
obtains  the  following  expression  for  the  displacement  response, 
in  terms  of  induced  displacement  input  d(t) : 


x(t)  =  e  b  "  (Acoscodt  +  Bsina)dt) 


+  C 


e  sin (a>dt  +  </>)  *d(t),  (11) 


where 


ft),  = 


=®„V i-o 

A  =  x(0  +), 


B  =  -^4r(0  +)+  — i(0  +)-^n-d(  0+), 


2$'®„ 


co. 


CO, 


C  =  — , 
cod 


CO, 


and 


)  =  tan 


1-V 


(12) 

(13) 

(14) 

(15) 

(16) 


The  asterisk  in  Equation  (11)  indicates  convolution,  defined  for 
arbitrary  functions  /  and  g  as  follows: 

f(t)*g(t):=  [  f(r)g(t-r)dr  .  (17) 

J0+ 

Convolution  term,  for  modal  (damped-sinusoid)  induced- 
disturbance  component  di  (t) 

To  proceed  it  will  be  necessary  to  evaluate  the  convolution 
term  of  Equation  (11)  for  an  exponentially  decaying  sinusoid 
dj(t)  [Eq.  (5)].  Expansion  of  the  integrand  and  application  of 
the  trigonometric  addition  formula 

sin(«  +  /?)  =  sinacos/?  +  cosasin/? ,  (18) 

yields 


C 


=  CDe  1  ( sii 


sini 


«.'/ .  /  .  ,\C'  ■ 


rW+$)J 


SI 


\  v  ‘  "JOf 
-CO  4/+^)Jt 

Integrating  [7];  and  making  the  substitutions 
M,  =  a,  -  qa)„  , 
Vi=<Od-(Oi, 
and  p,  =  a>d  +  ff); ; 

one  obtains  the  following  algebraic  result: 


n (a>dT+fa)  cos  cOjidr 


ea‘  sir{(OdT+t/))smcOjTdT) .  (19) 


(20) 

(21) 

(22) 


C 


e  b‘°"‘  s,m{a>dt  +  <j>)  *d[(t ) 


=CE|e_4  sin(ft}H-4) 


4  [/4-sin(v/+^)-v;cos(v;f+^)] 


2  2 
Bi  +v. 


4  4  sin(  Af +4~  Pi  coiQt + 0] 


-cos  {(Oft  +  fa) 


/H  sii^-q  cos^ 

A  sin^-p  cos^ 

> 

^4 +4  ] 

2(4  +4 ]  J 

eH‘  k  cos  (v;f  ■ +  4 +  V;  sin  + 4 

2 

(4 +4] 

4  k  cos  4 

I  +  !^) + A  sin  (a?  +  0] 

2|  4  +A 


A-cos^+v-sin^  /^-cos^+yCism^ 


2  2 
A  +L 


{4  +4  ] 


I  • (23) 


Define  now  the  following  variables,  for  convenience: 

Mi 


and 


K,;  = 


K 1,  = 


Kd;  =  ' 


2(h,2+p;)' 

V; 

2(p?+v2)' 

Pi' 

•(Mi2  +P2  ) 


(24) 

(25) 

(26) 
(27) 


Upon  substituting  from  Equations  (24)  through  (27),  Equation 
(23)  can  be  simplified  to  the  following  form: 

Ce  bW"r  sin (codt  +  4*  dj[t) 

=  C  Dte  'm"'  [  -K\j  cos  ( cod  t+</>  +  fa )+  x2i  cos  (codt  +  <j>-  fa ) 


-  kv  sin1 


in(cq/ 1  +  (j)  +  fa )  +  /r4lsin(fi)d  t  -fa )  ] 


-  C  Dte  a''  [  -K\i  cos  (cq  t  +  <j)i  +  (/> )+  K2j  cos  (cOjt  +  <f \  -  fa) 

-  K2i  sin  (cOjt  +  fa  +</>)- Kr4/sin(<y  /  +</> ;  -fa)  ] .  (28) 

Recall  [8]  the  trigonometric  relationship: 


^4cos9  +  RsinB  =  -^A2  +  B2  sin  0  +  tan  1 —  .  (29) 

B  , 
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With  the  help  of  Equation  (29),  Equation  (28)  can  now  be 
expressed  (after  some  algebraic  simplification),  as  the 
following  harmonic  sum: 


Ce  gc°"  sin(codt +  </>)*  dt(i) 

=  CDte  hC0"' [- 8V sin {wdt+<j)+<j)i+  9Xi ) 

+  £>2;Sin  (  G)dt  +  0  -  $  +  02i)  ] 

-CD/e  a,‘  [-  Su  sin  (cott  +(/>,  +  (j>  +  9Xi ) 

+  4, sin  ( (Ojt  +  </>i-(/  +  6^,)  ] ; 

where 

j~2  T  1 

2  yjMi+Vi 

2  2  1 

K2i  +  kM  = - 1 

\2  2 

2 \Mi  +Pi 


9U  =  tan  *  — -  =  tan 

Mi_ 

\vi. 

r  \ 

-  -1  K2i  -1 

92i  =  tan  -=*-  =  tan 

Ml 

km 

and 


9M  =  tan 


-1  K~, 


■  =  tan 


Mi 

~Pi 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


for  C  ,  p; ,  v,- ,  and  p;-  as  defined  by  Equations  (15),  (20),  (21), 
and  (22),  respectively. 


Position  response  x(f)  to  modally  componentiated  d(t ) 

Substitution  from  Equation  (4)  into  Equation  (11), 
followed  by  application  of  Equation  (30),  yields  the  following 
equation  for  the  SDOF  MSD  position  x(t) : 

x(t)=e  ?C°n  (Acos codt  +  B sin codt) 

+CX'D'i  sin(  t+<t>+<k  +6L) 

/= 1 

+S2isin{codt+<f-</i  +4,)] 

+e  "''[4sin(  (oit+(/)i+<p+9x) 

-  S2i  sinfrqt  +</>i~</>+  9M )]  | .  (36) 

Equation  (36)  describes  the  response  x{t)  to  the  total 
(composite)  induced  displacement  given  by  Equations  (4)  and 
(5). 

Relative  displacement  §(t),for  modally  componentiated  d(t) 

Substitution  from  Equations  (4),  (5),  and  (36)  into 
Equation  (6),  yields  the  following  expression  for  relative 
displacement,  of  the  SDOF  MSD  system  mass: 


8(t)  =  e  ?<°"  (Acosa>dt  +  Bsma>dt) 

+C±d\  e^'l-Sy  sin{o)dt+<j>+<l>i  +9y) 

i= 1 

+^,sin(©/+^-^  +8^j\ 

+e  a,[  4,sin(  wit+</>i  +</>+9x) 

- 4/  sixio^t  +  (j>i-(j>+ 93j ) 

-sin (o$f  +  $)/c]}.  (37) 

Velocity  response  x(t)  to  modally  componentiated  d{t) 

Differentiation  of  Equation  (36)  yields  the  following 
expression  for  the  velocity  response  i(t)  to  the  induced 
displacement  d(t): 

x (t) =-gaf,e  ^  (^4cos  codt  +Bsmcodt)  +  code~ia"'  (B cos  codt-A  sin  codt) 

V  [ 

+CYP\  -^^[-4- sin ((Qt-rf+fi  +6j,)  +<%/Sin(  codt +</>-$  +4)] 

i=l 

+o)deh°*[-8licos(cqt+</>+(/>,  +eti)+52icos(a)dt+<f>-fi  +4)] 

-ate  af[§dsir{coit+(/i  +</+0y)  -S2isir(o)it+^i  -fZ>+4,)] 

+a>j  e~°^[SyCos( wit-nf>j+ifH-9y)-S2icos( cOjt+tj)- 4-4)] .  (38) 

Acceleration  response  x(t)  to  modally  componentiated  d(t) 
Differentiating  Equation  (38)  with  respect  to  time,  and 
combining  terms,  yields  an  equation  for  the  acceleration 
response  x(t) : 

ix{t)=  (on-m^Y  (Acoscodt +Bsmcodt) 

-2 gconcode  ?  "  (B cos codt -Asmcodi] 

+cYPi \  (/ 4-4)e^[-4  si^®d^4+4')+4sin(®</+^+^)] 

+2 g«bO)de  ^[4-cos(cj  t+^+fa  +0li)-S2icos[(odt+</>-</>i  +4)] 

+(V  a/[4-sin(  cqt+tfi  +<p+d^j  -4jSir(  cqt+tfr  -</+92)\ 

+2 cqafi  a'' [— 4 cos {o)j  t+(/>/  +^+0u)+S2icos{ait+^  -fzl+4)]}.  (39) 

Note  that  the  number  of  terms  in  Equation  (39)  could  be 
reduced  by  a  factor  of  two,  for  computer  implementation,  by 
use  of  Equation  (29). 

Shock  response  spectra 

Now  it  is  possible  to  express  the  various  desired  shock 
response  spectra,  in  analytical  form.  The  spectral  displacement 
and  spectral  velocity  are  defined,  respectively,  by  Equations  (1) 
and  (2),  where  the  relative  displacement  5(t)  is  given  by 
Equation  (37).  The  spectral  acceleration  is  defined  by  Equation 
(3),  where  the  absolute  acceleration  x(t)  is  given  by  Equation 
(39). 
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SOLUTION  FOR  THE  UNDAMPED  CASE  ( q  =  0  ) 

Position  response  x(t) ,  undamped  case 

For  the  special  case  where  q  =  0 ,  the  position  response 
given  by  Equation  (36)  reduces  to  the  following: 


x(t)  =  Acosa>nt  +  B  sin  cont 

V 

(^,sin(<y,t+(4-  + 8 k) 
+S7i‘i\n(o)nt-<l>j  +4) 


+<a, 


!=1 


+e  4-sin  (  ®,++f4+4) 

-4  sin  (aft + fa +8^}}  ; 
where  4/  =  — .  , 

2  V«/  +(®«-®/)2 
1 

°2i  -  I  2  ’ 

2tJ  a,  +  (con+coj)2 

/  \ 

-3  +  ai 

8U  =  tan  - - —  , 

-if  ai 

4,-  =  tan  - - —  , 

and  8m  =  tan  '  f - — - | . 


a, 


a, 


v®«  +®.v 


a, 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


Relative  displacement  8(f),  undamped  case 

With  no  damping,  the  relative  displacement  given  by 
Equation  (37)  reduces  to  the  following: 

(?(t)=  A  cos  a> nt  +  Borneo nt 


V 

+C0tfYJD,{  -4- sin t+$  +4)  +4/sin(®„t-$  +80J 
1=1 

+e  [  4-sin  +4)  -  4 sin(qf +$  +4,) 

-sin  (®,f  +  (4  )/»«]}>  (46) 


where  the  constants  4  and  4  can  6c  appropriately  simplified, 
as  given  in  Equations  (41)  through  (45). 

Acceleration  response  x(f),  undamped  case 

Without  damping,  the  absolute  acceleration  given  by 
Equation  (39)  reduces  (again,  as  above,  with  simplified 
constants)  to  the  following  form: 


*W=  -4  (Acoscont +Bsmcont) 

+®„  [-4 sin {co„t+<k  +4  )+4sin(f(j,/-^  +4)] 

i= 1 

+{ai  a,'[Sysir{(oit+(/>i  +4)  -4;sin(®,.t+^  +4,.)] 

+ 2a; coje^,> [- 4 cos ( t+4  +4)+4cos((w,-t+4  +4,)]}.  (47) 

VERIFICATION  AND  MATLAB  IMPLEMENTATION 

An  algebraic  check  of  Equations  (37)  and  (39) 

In  the  limiting  case,  where  g  =  0 ,  Equation  (8)  can  be 
expressed  by 

x  +  ro'S  =  0 ,  (48) 

where  8  is  defined  by  Equation  (6).  The  expressions  for  8 
and  X  from  Equations  (37)  and  (39),  respectively,  should 
satisfy  Equation  (48)  identically.  Substituting,  one  obtains 

x+co~S=  -co2„(Acoscont+Bsmcont) 

V  ( 

4  [4  sin [cont+<k  +4  )-4sin(ffl,/-$  +4)] 

i=i 

+{4  -cf  "''[4-sin (ajt+fi  +4) -S2isir((Oit+</>i  +4,)] 

+2 aicoie~ai‘  [-4  cos(  cq  t+tfr  +4)+4  cos(  °¥ +<fi  +4)  ]} 

+  co~(Acosa>„t  +  Bsmcont) 

+°4YjDi<°n{  “4  sin(^  ?+(4  +61/)  +4-sin(®„t-4  +4) 

1=1 

+e^'[4-sin(£qt+(4  +4)  -4-sin(fi^+$+4) 


-sin (qf  +  $)/®„]}, 


(49) 


which  reduces  readily  to 

2  „ 

x +  £»,,£>  = 

(On  [4  -4  +4  V^[4sin(  cqt-uj  +4)  -4,sin(  cqt^  +4) ] 

-2 OtCDf*  [4  cos  (q  t+</>i  +8y) -4  cos  ( cof  +<t>i  +4)  ] 
-®„e~^sin(qt+f4)}.  (50) 
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Further  expansion  of  terms  yields  the  following: 
x  +  ro„o  = 


-d2i  cos03isin  [cott +</>,)- d2i  sin03i  cosRr+$)] 


-  2 aieh  R  cos^  cos  Rf  +  fai )  -  Su  sin#1(  sin  Rf  +  fat ) 

-  52i  cos#3,'  cos Rf +</>i)  +  52j  sin#,,  sin  R/  +  faj )] 

-co„  sinRf  +  +/)}■  (51) 

Combination  of  like  trigonometric  terms  yields 

2f 

X  +  K)„0  = 

A  e~a,t {sinRf  +  fa)  [R  cos^,  -S2i cos 62l){a)n  -cof  +aj  j 

1=1 

+R  sin^i,  -4  sin6>3,.)(2«i.^.)-fflj 

+  cosRf  +  fa )  [R  sin^  -  4  sin4  ){con  -  a>-  +  a]  j 

-Rcos4  ~4  cos4)(2a,. «,.)]}.  (52) 

In  order  for  Equation  (48)  to  hold  nontrivially,  it  is  necessary 
that  the  square -bracketed  terms  in  Equation  (52)  be  zero,  for  all 
indices  i.  In  particular,  it  is  necessary  that 

[R  cos4  -4  cos03i)(al  -a-  +a-  j 

+R  sin4  -4,  si  n  6*,,)  {la[  a>i  )-«„]  =  0 ,  (53) 

and  that 

[R  sin^i,  “4  sin4)^  ~0)f  +aj  j 

-Rcos4  - 4  co sd2j ) (2a,- &>,)]}  =  0 .  (54) 

From  Equations  (43)  and  (45),  one  can  obtain  the  following: 


and 


sin  4  = 


cos  4  = 


sin  4;  = 


cos  4-  = 


(55) 

(56) 

(57) 

(58) 


Substitution  from  Equations  (41),  (42),  and  (55)  through  (58), 
into  Equations  (53)  and  (54),  proves  the  latter  two  equations  to 
express  valid  identities.  This  in  turn  demonstrates  that 
Equations  (46)  and  (47)  satisfy  Equation  (48)  identically 
(Q.E.D.). 


A  numerical  check  of  Equation  (37) 

To  verify  the  results  for  the  more  general,  underdamped 
case,  the  algebraic  equation  for  relative  displacement  5(f) 
[Equation  (37)]  was  first  implemented  in  MATLAB  code. 
Then,  for  selected  damping  ratios  (i.e.,  of  the  hypothetical 
SDOF  MSD  system),  and  for  various  induced-displacement 
inputs  [Eqs.  (4)  and  (5)],  numerical  evaluations  were  made  of 
the  spectral  displacement,  using  Equation  (1)  with  the  algebraic 
results  given  by  Equation  (37).  These  results  were  compared 
with  numerical  evaluations  of  the  same  spectral  quantity  (using 
MATLAB)  with  the  convolution  integrals  determined  by  direct 
numerical  integration  [Eqs.  (1),  (4),  (5),  and  (11)].  The  results 
from  algebraic  substitution  and  numerical  integration  are 
identical. 

For  a  representative  example,  consider  an  induced- 
displacement  input  with  components  described  by  the 
parameters  of  Table  1.  Table  2  displays  the  results  of 
calculations  of  corresponding  spectral  displacement  by  the  two 
methods  described  above.  The  two  approaches  are  seen  to  give 
identical  results. 

Table  1.  Parameters  for  Induced-Displacement 
Input  Components  dt(t) 


i  =  1 

i  =  2 

i  =  3 

i  =  4 

i  =  5 

A 

0.009 

0.002 

0.002 

0.001 

0.006 

a  ,■ 

0.50 

0.25 

0.15 

0.95 

0.75 

CD, 

250 

100 

750 

900 

500 

<t>i 

0 

0 

0 

0 

0 

Table  2.  Comparison  of  Spectral  Displacements 
Evaluated  by  Different  Techniques  (time  step  0.001  s) 


SRS  Freq.  oj„ 

(<7  =  0.1) 

AlRi)  from 
Algebraic  Substitution 

4(®h)  from 
Numerical  Integration 

1 

0.01639525625157 

0.01639525625157 

2 

0.01641587042504 

0.01641587042504 

3 

0.01644317777833 

0.01644317777833 

4 

0.01647694228270 

0.01647694228270 

5 

0.01651690741770 

0.01651690741770 

6 

0.01656279716325 

0.01656279716325 

7 

0.01661431703708 

0.01661431703708 

8 

0.01667115517439 

0.01667115517439 

9 

0.01673298344652 

0.01673298344652 

10 

0.01679945861516 

0.01679945861516 
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Figure  3  presents  plots  of  the  spectral  quantities,  for  this  input 
d(t)  (Table  1),  and  for  this  damping  ratio  g  (Table  2),  with 
zero  initial  conditions. 


Figure  3.  Shock-Response  Spectra  Using  the  Induced- 
Displacement  Input  of  Table  1,  for  damping  ratio  0.1. 


CONCLUSION 

This  paper  has  presented  analytical  equations  describing 
the  spectral  displacement  (displacement  SRS),  the  spectral 
velocity  (pseudovelocity  SRS),  and  the  spectral  acceleration 
(absolute  acceleration  SRS),  for  a  proportionally  damped,  linear 
system.  For  such  systems  an  impulsive  mechanical-shock 
disturbance  produces  a  vibratory  response  expressible 
analytically  in  terms  of  the  system  modes  and  modeshapes, 
using  normal-mode  orthogonalization.  In  particular,  the 
response  has  the  form  of  a  linear  combination  of  exponentially 
decaying  sinusoids,  of  various  amplitudes  and  phase  shifts.  A 
response  of  such  a  form  can  be  represented  by  SRS’s  for  which 
this  paper  has  provided  analytical  descriptions.  The 
displacement  and  pseudovelocity  SRS’s  are  defined, 
respectively,  by  Equations  (1)  and  (2),  where  the  relative 
displacement  8{t)  is  given  analytically  by  Equation  (37).  The 
spectral  acceleration  is  defined  by  Equation  (3),  where  the 
absolute  acceleration  is  given  analytically  by  Equation  (39). 

Having  these  analytical  expressions  for  the  various  SRS’s 
permits  the  SRS’s  to  be  computed  exactly  for  impulsive  shock 
disturbances,  without  necessitating  numerical  evaluation  of  the 
convolution  integral.  The  analytical  expressions  can  be  used 
for  other,  non-impulsive  input  loads,  even  those  for  which  there 
is  no  simple  analytical  description,  provided  they  can  be 


approximated  as  impulses.  The  physical  system  itself  serves  as 
a  modal  filter  of  the  shock  input,  to  produce  a  vibratory 
response,  known  at  any  desired  point  on  the  system  in  terms  of 
its  exponentially  decaying  sinusoidal  components.  Each  of 
these  components  is  known  in  terms  of  four  values:  its 
frequency,  its  decay  rate  (time  constant),  its  phase  angle,  and  its 
amplitude.  From  these  values  the  desired  SRS  can  be 
determined  by  evaluating,  at  each  point  of  a  discretized 
continuum  of  frequencies  (i.e.,  those  of  the  conceptualized 
SDOF  MSD  systems),  the  maximum  (or  minimum)  of  a  time 
function  consisting  of  simple  algebraic  expressions  involving 
simple  trigonometric  operations.  No  numerical  convolution  is 
required,  because  the  integrations  have  already  been 
accomplished  analytically  (and  exactly,  for  impulse  loading), 
with  the  results  incorporated  into  the  algebraic  expressions. 
This  method  can  provide  for  accurate  SRS  computation 
irrespective  of  the  input  shock’s  exact  shape,  provided  the  input 
is  approximately  impulsive.  For  linear  systems  the  described 
method  can  be  used  as  a  benchmark  to  evaluate  the  accuracy  of 
other  methods  of  SRS  determination.  It  can  also  be  used  to 
determine  the  minimum  number  of  modes  required,  in  a 
system’s  finite-element  model,  to  produce  an  SRS  of  specified 
accuracy. 
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US  ARMY  ARDEC 
AMSTA  AR  CCH  B 
P  DONADIA 
F  DONLON 
P  VALENTI 
C  KNUTSON 
G  EUSTICE 
K  HENRY 
J  MCNABOC 
G  WAGNECZ 
R  SAYER 
F  CHANG 

PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  CCL 
F  PUZYCKI 
R  MCHUGH 
D  CONWAY 
E  JAROSZEWSKI 
R  SCHLENNER 
M  CLUNE 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  PM  ARMS 

SFAE  GCSS  ARMS 
BLDG  171 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  WEA 
J  BRESCIA 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  PM  MAS 

SFAE  AMO  MAS 
PICATINNY  ARSENAL  NJ 
07806-5000 
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NO.  OF 

COPIES  ORGANIZATION 
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COPIES  ORGANIZATION 


1  PM  MAS 

SFAE  AMO  MAS 
CHIEF  ENGINEER 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  PM  MAS 

SFAE  AMO  MAS  PS 
PICATINNY  ARSENAL  NJ 
07806-5000 

2  PM  MAS 

SFAE  AMO  MAS  LC 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
PRODUCTION  BASE 
MODERN  ACTY 
AMSMC  PBM  K 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  TACOM 
PM  COMBAT  SYSTEMS 
SFAE  GCS  CS 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

US  ARMY  TACOM 
AMSTA  SF 

WARREN  MI  48397-5000 

1  DIRECTOR 

AIR  FORCE  RESEARCH  LAB 
MLLMD 
D  MIRACLE 
2230  TENTH  ST 

WRIGHT  PATTERSON  AFB  OH 
45433-7817 

1  OFC  OF  NAVAL  RESEARCH 
J  CHRISTODOULOU 
ONR  CODE  332 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5600 

1  US  ARMY  CERL 
R  LAMPO 

2902  NEWMARK  DR 
CHAMPAIGN  IL  61822 


1  COMMANDER 

US  ARMY  TACOM 
PM  SURVIVABLE  SYSTEMS 
SFAE  GCSS  W  GSI H 
M  RYZYI 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

US  ARMY  TACOM 
CHIEF  ABRAMS  TESTING 
SFAE  GCSS  W  AB  QT 
T  KRASKIEWICZ 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

WATERVLIET  ARSENAL 
SMCWV  QAE  Q 
B  VANINA 
BLDG  44 

WATERVLIET  NY  12189-4050 

1  TNG,  DOC,  &  CBT  DEV 
ATZK  TDD  IRSA 

A  POMEY 
FT  KNOX  KY  40121 

2  HQ  IOC  TANK 
AMMUNITION  TEAM 
AMSIO  SMT 

R  CRAWFORD 
W  HARRIS 

ROCK  ISLAND  IL  61299-6000 

2  COMMANDER 

US  ARMY  AMCOM 
AVIATION  APPLIED  TECH  DIR 
J  SCHUCK 

FT  EUSTIS  VA  23604-5577 

1  NSWC 

DAHLGREN  DIV  CODE  G06 
DAHLGREN  VA  22448 

2  US  ARMY  CORPS  OF  ENGR 
CERDC 

T  LIU 
CEW  ET 
T  TAN 

20  MASSACHUSETTS  AVE  NW 
WASHINGTON  DC  20314 
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COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  US  ARMY  COLD  REGIONS 
RSCH  &  ENGRNG  LAB 
P  DUTTA 
72  LYME  RD 
HANOVER  NH  03755 

14  COMMANDER 

US  ARMY  TACOM 
AMSTA  TR  R 
R  MCCLELLAND 
D  THOMAS 
J  BENNETT 
D  HANSEN 
AMSTA  JSK 
S  GOODMAN 
J  FLORENCE 
K  IYER 

D  TEMPLETON 
A  SCHUMACHER 
AMSTA  TR  D 
D  OSTBERG 
L  HINOJOSA 
B  RAJU 
AMSTA  CS  SF 
H  HUTCHINSON 
F  SCHWARZ 
WARREN  MI  48397-5000 

14  BENETLABS 
AMSTA  AR  CCB 
R  FISCELLA 
M  SOJA 
E  KATHE 
M  SCAVULO 
G  SPENCER 
P  WHEELER 
S  KRUPSKI 
J  VASILAKIS 
G  FRIAR 
R  HASENBEIN 
AMSTA  CCB  R 
S  SOPOK 
E  HYLAND 
D  CRAYON 
R  DILLON 

WATERVLIETNY  12189-4050 

1  USA  SBCCOM  PM  SOLDIER  SPT 
AMSSB  PM  RSS  A 
J  CONNORS 
KANSAS  ST 
NATICK  MA  01760-5057 


1  NSWC 

TECH  LIBRARY  CODE  323 
17320  DAHLGRENRD 
DAHLGREN  VA  22448 

2  USA  SBCCOM 
MATERIAL  SCIENCE  TEAM 
AMSSB  RSS 

J  HERBERT 
M  SENNETT 
KANSAS  ST 
NATICK  MA  01760-5057 

2  OFC  OF  NAVAL  RESEARCH 
D  SIEGEL  CODE  351 
J  KELLY 

800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 

1  NSWC 

CRANE  DIVISION 
M  JOHNSON  CODE  20H4 
LOUISVILLE  KY  40214-5245 

2  NSWC 

U  SORATHIA 
C  WILLIAMS  CD  6551 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817 

2  COMMANDER 
NSWC 

CARDEROCK  DIVISION 
R  PETERSON  CODE  2020 
M  CRITCHFIELD  CODE  1730 
BETHESDA  MD  20084 

8  DIRECTOR 

US  ARMY  NGIC 
D  LEITER  MS  404 
MHOLTUS  MS  301 
M  WOLFE  MS  307 
S  MINGLEDORF  MS  504 
J  GASTON  MS  301 
W  GSTATTENBAUER  MS  304 
R  WARNER  MS  305 
J  CRIDER  MS  306 
2055  BOULDERS  RD 
CHARLOTTESVILLE  VA 
22911-8318 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  NAVAL  SEA  SYSTEMS  CMD 
D  LIESE 

1333  ISAAC  HULL  AVE  SE  1100 
WASHINGTON  DC  20376-1100 

1  EXPEDITIONARY  WARFARE 
DIV  N85 
F  SHOUP 

2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

8  US  ARMY  SBCCOM 

SOLDIER  SYSTEMS  CENTER 
BALLISTICS  TEAM 
J  WARD 
WZUKAS 
P  CUNNIFF 
J  SONG 

MARINE  CORPS  TEAM 

J  MACKIEWICZ 

BUS  AREA  ADVOCACY  TEAM 

W  HASKELL 

AMSSB  RCP  SS 

W  NYKVIST 

S  BEAUDOIN 

KANSAS  ST 

NATICK  MA  01760-5019 

7  US  ARMY  RESEARCH  OFC 
A  CROWSON 
H  EVERETT 
J  PRATER 
G  ANDERSON 
D  STEPP 
D  KISEROW 
J  CHANG 
POBOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

1  AFRL  MLBC 

2941  PSTRM  136 

WRIGHT  PATTERSON  AFB  OH 

45433-7750 

1  DIRECTOR 

LOS  ALAMOS  NATL  LAB 
F  L  ADDESSIO  T  3  MS  5000 
PO  BOX  1633 
LOS  ALAMOS  NM  87545 


8  NSWC 

J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
R  D  COOPER  CODE  G32 
J  FRAYSSE  CODE  G33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGREN  VA  22448 

1  NSWC 

CARDEROCK  DIVISION 
R  CRANE  CODE  6553 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817-5700 

1  AFRL  MLSS 
R  THOMSON 

2179  12TH  ST  RM  122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 

2  AFRL 

F  ABRAMS 

J  BROWN 

BLDG  653 

2977  P  ST  STE  6 

WRIGHT  PATTERSON  AFB  OH 

45433-7739 

5  DIRECTOR 
LLNL 

R  CHRISTENSEN 
S  DETERESA 
F  MAGNESS 
M  FINGER  MS  313 
M  MURPHY  L  282 
PO  BOX  808 
LIVERMORE  CA  94550 

1  AFRL  MLS  OL 
L  COULTER 
5851  F  AVE 
BLDG  849  RM  ADI  A 
HILL  AFB  UT  84056-5713 

1  OSD 

JOINT  CCD  TEST  FORCE 
OSD  JCCD 
R  WILLIAMS 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  29180-6199 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


3  DARPA 

M  VANFOSSEN 
S  WAX 

L  CHRISTODOULOU 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

2  SERDP  PROGRAM  OFC 
PM  P2 

C  PELLERIN 
B  SMITH 

901  N  STUART  ST  STE  303 
ARLINGTON  VA  22203 

1  OAK  RIDGE  NATL  LAB 
R  M  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  3  783 1  -6 1 95 

1  OAK  RIDGE  NATL  LAB 
C  EBERLE  MS  8048 
PO  BOX  2008 
OAK  RIDGE  TN  3 7831 

3  DIRECTOR 
SANDIA  NATL  LABS 
APPLIED  MECHS  DEPT 
MS  9042 

J  HANDROCK 
Y  R  KAN 
J  LAUFFER 
PO  BOX  969 

LIVERMORE  CA  94551-0969 

1  OAK  RIDGE  NATL  LAB 
C  D  WARREN  MS  8039 
PO  BOX  2008 
OAK  RIDGE  TN  37831 

4  NIST 

M  VANLANDINGHAM  MS  8621 
J  CHIN  MS  8621 
J  MARTIN  MS  8621 
D  DUTHINH  MS  8611 
100  BUREAU  DR 
GAITHERSBURG  MD  20899 

1  HYDROGEOLOGIC  INC 
SERDP  ESTCP  SPT  OFC 
S  WALSH 

1155  HERNDON  PKWY  STE  900 
HERNDON  VA  20170 


3  NASA  LANGLEY  RESEARCH  CTR 
AMSRD  ARL  VS 
W  ELBER  MS  266 
F  BARTLETT  JR  MS  266 
G  FARLEY  MS  266 
HAMPTON  VA  23681-0001 

1  NASA  LANGLEY  RESEARCH  CTR 
T  GATES  MS  188E 
HAMPTON  VA  23661-3400 

1  FHWA 

E  MUNLEY 

6300  GEORGETOWN  PIKE 
MCLEAN  VA  22101 

1  USDOT  FEDERAL  RAILROAD 
M  FATEH  RDV  31 
WASHINGTON  DC  20590 

3  CYTEC  FIBERITE 
R DUNNE 
D  KOHLI 
R  MAYHEW 
1300  REVOLUTION  ST 
HAVRE  DE  GRACE  MD  21078 

1  DIRECTOR 
NGIC 
IANG  TMT 
2055  BOULDERS  RD 
CHARLOTTESVILLE  VA 
22911-8318 

1  SIOUX  MFG 
B  KRIEL 
PO  BOX  400 

FT  TOTTEN  ND  58335 

2  3  TEX  CORP 

A  BOGDANOVICH 
J  SINGLETARY 
109  MACKENAN  DR 
CARY  NC  27511 

1  3M  CORP 
J  SKILDUM 

3M  CENTER  BLDG  60  IN  01 
ST  PAUL  MN  55144-1000 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  DIRECTOR 

DEFENSE  INTLLGNC  AGNCY 
TA  5 

K  CRELLING 
WASHINGTON  DC  20310 

1  ADVANCED  GLASS  FIBER  YARNS 
T  COLLINS 

281  SPRING  RUN  LANE  STE  A 
DOWNINGTON  PA  19335 

1  COMPOSITE  MATERIALS  INC 
D  SHORTT 
19105  63  AVE  NE 
PO  BOX  25 

ARLINGTON  WA  98223 

1  JPS  GLASS 
L  CARTER 
PO  BOX  260 
SLATER  RD 
SLATER  SC  29683 

1  COMPOSITE  MATERIALS  INC 
R  HOLLAND 
1 1  JEWEL  CT 
ORINDA  CA  94563 

1  COMPOSITE  MATERIALS  INC 
C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 

2  SIMULA 

J  COLTMAN 
R  HUYETT 
10016  S51ST  ST 
PHOENIX  AZ  85044 

2  PROTECTION  MATERIALS  INC 
M  MILLER 
F  CRILLEY 
14000  NW  58  CT 
MIAMI  LAKES  FL  33014 

2  FOSTER  MILLER 
M  ROYLANCE 
WZUKAS 
195  BEAR  HILL  RD 
WALTHAM  MA  02354-1196 


1  ROM  DEVELOPMENT  CORP 
R  O  MEARA 

136  SWINEBURNE  ROW 
BRICK  MARKET  PLACE 
NEWPORT  RI  02840 

2  TEXTRON  SYSTEMS 
T  FOLTZ 

M  TREASURE 
1449  MIDDLESEX  ST 
LOWELL  MA  01851 

1  O  GARA  HESS  &  EISENHARDT 
M  GILLESPIE 

9113  LESAINT  DR 
FAIRFIELD  OH  45014 

2  MILLIKEN  RESEARCH  CORP 
H  KUHN 

M  MACLEOD 
PO  BOX  1926 
SPARTANBURG  SC  29303 

1  CONNEAUGHT  INDUSTRIES  INC 
J  SANTOS 
PO  BOX  1425 
COVENTRY  RI  02816 

1  ARMTEC  DEFENSE  PRODUCTS 
S  DYER 
85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 

1  NATL  COMPOSITE  CTR 
T  CORDELL 
2000  COMPOSITE  DR 
KETTERING  OH  45420 

3  PACIFIC  NORTHWEST  LAB 
M  SMITH 

G  VAN  ARSDALE 
R  SHIPPELL 
PO  BOX  999 
RICHLAND  WA  99352 

1  SAIC 

M  PALMER 

1410  SPRING  HILL  RD  STE  400 
MS  SH4 5 

MCLEAN  VA  22 102 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  ALLIANT  TECHSYSTEMS  INC 

4700  NATHAN  LN  N 
PLYMOUTH  MN  55442-2512 

1  APPLIED  COMPOSITES 
W  GRISCH 

333  NORTH  SIXTH  ST 
ST  CHARLES  IL  60174 

1  CUSTOM  ANALYTICAL 
ENG  SYS  INC 
A  ALEXANDER 
13000  TENSOR  LANE  NE 
FLINTSTONE  MD  21530 

1  AAI CORP 

DR  N  B  MCNELLIS 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

1  OFC  DEPUTY  UNDER  SEC  DEFNS 
J  THOMPSON 

1745  JEFFERSON  DAVIS  HWY 
CRYSTAL  SQ4STE501 
ARLINGTON  VA  22202 

3  ALLIANT  TECHSYSTEMS  INC 

J  CONDON 
E  LYNAM 
J  GERHARD 
WV01  16  STATE  RT  956 
PO  BOX  210 
ROCKET  CENTER  WV 
26726-0210 

1  PROJECTILE  TECHNOLOGY  INC 
515  GILES  ST 

HAVRE  DE  GRACE  MD  21078 

1  HEXCEL  INC 
RBOE 

POBOX  18748 

SALT  LAKE  CITY  UT  841 18 

1  PRATT  &  WHITNEY 
C  WATSON 

400  MAIN  ST  MS  114  37 
EAST  HARTFORD  CT  06108 


5  NORTHROP  GRUMMAN 
B  IRWIN 
K  EVANS 
D  EWART 

A  SHREKENHAMER 
J  MCGLYNN 
BLDG  160  DEPT  3700 
1 100  WEST  HOLLYVALE  ST 
AZUSA  CA  91701 

1  HERCULES  INC 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  BRIGS  COMPANY 
J  BACKOFEN 
2668  PETERBOROUGH  ST 
HERNDON  VA  22071-2443 

1  ZERNOW  TECHNICAL  SERVICES 
L  ZERNOW 

425  W  BONITA  AVE  STE  208 
SAN  DIMAS  CA  91773 

1  GENERAL  DYNAMICS  OTS 
L  WHITMORE 

10101  NINTH  ST  NORTH 
ST  PETERSBURG  FL  33702 

2  GENERAL  DYNAMICS  OTS 
FLINCHBAUGH  DIV 

K  LINDE 
T  LYNCH 
PO  BOX  127 
RED  LION  PA  17356 

1  GKN  WESTLAND  AEROSPACE 
D  OLDS 

450  MURDOCK  AVE 
MERIDEN  CT  06450-8324 

2  BOEING  ROTORCRAFT 
P  MINGURT 

P  HANDEL 

800  B  PUTNAM  BLVD 

WALLINGFORD  PA  19086 
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COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


5  SIKORSKY  AIRCRAFT 
G  JACARUSO 
T  CARSTENSAN 
B  KAY 

S  GARBO  MS  S330A 
J  ADELMANN 
6900  MAIN  ST 
PO  BOX  9729 

STRATFORD  CT  06497-9729 

1  AEROSPACE  CORP 
G  HAWKINS  M4  945 

2350  E  EL  SEGUNDO  BLVD 
EL  SEGUNDO  CA  90245 

2  CYTEC  FIBERITE 
M  LIN 

W  WEB 

1440  N  KRAEMER  BLVD 
ANAHEIM  CA  92806 

2  UDLP 

G  THOMAS 
M  MACLEAN 
PO  BOX  58123 
SANTA  CLARA  CA  95052 

1  UDLP  WARREN  OFC 
ALEE 

31201  CHICAGO  RD  SOUTH 
SUITE  B 102 
WARREN  MI  48093 

2  UDLP 

R  BRYNSVOLD 
P  JANKEMS  170 
4800  EAST  RIVER  RD 
MINNEAPOLIS  MN  55421-1498 

1  LOCKHEED  MARTIN 
SKUNK  WORKS 
D  FORTNEY 
1011  LOCKHEED  WAY 
PALMDALE  CA  93599-2502 

1  LOCKHEED  MARTIN 
R  FIELDS 
5537  PGA  BLVD 
SUITE  4516 
ORLANDO  FL  32839 


1  NORTHRUP  GRUMMAN  CORP 
ELECTRONIC  SENSORS 
&  SYSTEMS  DIV 
E  SCHOCH  MS  V  16 
1745 AW  NURSERY  RD 
LINTHICUM  MD  21090 

1  GDLS  DIVISION 
D  BARTLE 
POBOX  1901 
WARREN  MI  48090 

2  GDLS 
DREES 
M  PASIK 
PO  BOX  2074 
WARREN  MI  48090-2074 

1  GDLS 

MUSKEGON  OPER 
M  SOIMAR 
76  GETTY  ST 
MUSKEGON  MI  49442 

1  GENERAL  DYNAMICS 
AMPHIBIOUS  SYS 
SURVIVABILITY  LEAD 
G  WALKER 
991  ANNAPOLIS  WAY 
WOODBRIDGE  VA  22191 

6  INST  FOR  ADVANCED 
TECH 
H  FAIR 
I  MCNAB 
P  SULLIVAN 
S  BLESS 
W  REINECKE 
C  PERSAD 

3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

1  ARROW  TECH  AS  SOC 

1233  SHELBURNE  RD  STE  D8 
SOUTH  BURLINGTON  VT 
05403-7700 

1  R  EICHELBERGER 
CONSULTANT 
409  W  CATHERINE  ST 
BEL  AIR  MD  21014-3613 
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2 

1 

1 

1 

1 

1 
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COPIES  ORGANIZATION 


SAIC 

G  CHRYSSOMALLIS 

8500  NORMANDALE  LAKE  BLVD 

SUITE  1610 

BLOOMINGTON  MN  55437-3828 

UCLA  MANE  DEPT  ENGR IV 
H  T  HAHN 

LOS  ANGELES  CA  90024-1597 

UNIV  OF  DAYTON 
RESEARCH  INST 
R  Y  KIM 
AKROY 

300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0168 

UMASS  LOWELL 
PLASTICS  DEPT 
N  SCHOTT 
1  UNIVERSITY  AVE 
LOWELL  MA  01854 

IIT  RESEARCH  CTR 
DROSE 
201  MILL  ST 
ROME  NY  13440-6916 

GA  TECH  RESEARCH  INST 
GA  INST  OF  TCHNLGY 
P  FRIEDERICH 
ATLANTA  GA  30392 

MICHIGAN  ST  UNIV 
MSM  DEPT 
R  AVERILL 
3515  EB 

EAST  LANSING  MI  48824-1226 

UNIV  OF  WYOMING 
D  ADAMS 
PO  BOX  3295 
LARAMIE  WY  82071 

PENN  STATE  UNIV 
R  S  ENGEL 

245  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  16801 


2  PENN  STATE  UNIV 
R  MCNITT 
C  BAKIS 

212  EARTH  ENGR 
SCIENCES  BLDG 
UNIVERSITY  PARK  PA  16802 

1  PURDUE  UNIV 

SCHOOL  OF  AERO  &  ASTRO 
C  T  SUN 

W  LAFAYETTE  IN  47907-1282 

1  STANFORD  UNIV 

DEPT  OF  AERONAUTICS 
&  AEROBALLISTICS 
S  TSAI 

DURANT  BLDG 
STANFORD  CA  94305 

1  UNIV  OF  MAINE 

ADV  STR  &  COMP  LAB 
R  LOPEZ  ANIDO 
5793  AEWC  BLDG 
ORONO  ME  04469-5793 

1  JOHNS  HOPKINS  UNIV 
APPLIED  PHYSICS  LAB 
P  WIENHOLD 
11100  JOHNS  HOPKINS  RD 
LAUREL  MD  20723-6099 

1  UNIV  OF  DAYTON 
J  M  WHITNEY 
COLLEGE  PARK  AVE 
DAYTON  OH  45469-0240 

1  NORTH  CAROLINA  ST  UNIV 
CIVIL  ENGINEERING  DEPT 
W  RASDORF 
PO  BOX  7908 
RALEIGH  NC  27696-7908 

5  UNIV  OF  DELAWARE 

CTR  FOR  COMPOSITE  MTRLS 
J  GILLESPIE 
M  SANT ARE 
S  YARLAGADDA 
S  ADVANI 
D  HEIDER 
201  SPENCER  LAB 
NEWARK  DE  19716 
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COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  DEPT  OF  MTRLS 
SCIENCE  &  ENGRG 
UNIV  OF  ILLINOIS 
AT  URBANA  CHAMPAIGN 
J  ECONOMY 

1304  WEST  GREEN  ST  115B 
URBANA  IL  61801 

1  UNIV  OF  MARYLAND 

DEPT  OF  AEROSPACE  ENGRG 
A  J  VIZZ1NI 

COLLEGE  PARK  MD  20742 

1  DREXEL  UNIV 

A  S  D  WANG 
3141  CHESTNUT  ST 
PHILADELPHIA  PA  19104 

3  UNIV  OF  TEXAS  AT  AUSTIN 
CTR  FOR  ELECTROMECHANICS 
J  PRICE 
A  WALLS 
J  KITZMILLER 
10100  BURNET  RD 
AUSTIN  TX  78758-4497 

3  VA  POLYTECHNICAL 

INST  &  STATE  UNIV 
DEPT  OF  ESM 
M  W  HYER 
K  REIFSNIDER 
R  JONES 

BLACKSBURG  VA  24061-0219 

1  SOUTHWEST  RESEARCH  INST 
ENGR  &  MATL  SCIENCES  DIV 
J  RIEGEL 

6220  CULEBRA  RD 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78228-0510 

1  BATELLE  NATICK  OPERS 
B  HALPIN 
313  SPEENST 
NATICK  MA  01760 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  WM  MB 
A  FRYDMAN 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


ABERDEEN  PROVING  GROUND 

1  US  ARMY  ATC 

CSTE  DTC  AT  AC  I 
W  C  FRAZER 
400  COLLERAN  RD 
APG  MD  21005-5059 

91  DIR  US  ARL 

AMSRD  ARL  Cl 
AMSRD  ARL  O  AP  EG 
M  ADAMSON 
AMSRD  ARL  SL  BA 
AMSRD  ARL  SL  BB 
D  BELY 

AMSRD  ARL  WM 
J  SMITH 
H  WALLACE 
AMSRD  ARL  WM  B 
A  HORST 
T  KOGLER 
AMSRD  ARL  WM  BA 
D  LYON 

AMSRD  ARL  WM  BC 
J  NEWILL 
PPLOSTINS 
A  ZIELINSKI 
AMSRD  ARL  WM  BD 
P  CONROY 
B  FORCH 
M  LEADORE 
C  LEVERITT 
RLIEB 

R  PESCE  RODRIGUEZ 
B  RICE 

AMSRD  ARL  WM  BF 
S  WILKERSON 
AMSRD  ARL  WM  M 
B  FINK 
J  MCCAULEY 
AMSRD  ARL  WM  MA 
L  GHIORSE 
S  MCKNIGHT 
E  WETZEL 
AMSRD  ARL  WM  MB 
J  BENDER 
T  BOGETTI 
L  BURTON 
R  CARTER 
KCHO 

W  DE  ROSSET 
G  DEWING 
R  DOWDING 
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W  DRYSDALE 

AMSRD  ARL  WM  TB 

R  EMERSON 

P  BAKER 

D  HENRY 

AMSRD  ARL  WM  TC 

D  HOPKINS 

R  COATES 

RKASTE 

AMSRD  ARL  WM  TD 

L  KECSKES 

D  DANDEKAR 

M  MINNICINO 

T  HADUCH 

B  POWERS 

T  MOYNIHAN 

D  SNOHA 

M  RAFTENBERG 

J  SOUTH 

S  SCHOENFELD 

M  STAKER 

T  WEERASOORIY A 

J  SWAB 

AMSRD  ARL  WM  TE 

J  TZENG 

A  NIILER 

AMSRD  ARL  WM  MC 

J  POWELL 

J  BEATTY 
R  BOSSOLI 
ECHIN 

S  CORNELISON 
D  GRANVILLE 
B  HART 
J  LASALVIA 
J  MONTGOMERY 
F  PIERCE 
ERIGAS 
W  SPURGEON 

AMSRD  ARL  WM  MD 
B  CHEESEMAN 
P  DEHMER 
R  DOOLEY 
G  GAZONAS 
S  GHIORSE 
C  HOPPEL 
M  KLUSEWITZ 
WROY 
J  SANDS 
D  SPAGNUOLO 
S  WALSH 
S  WOLF 

AMSRD  ARL  WM  RP 
J  BORNSTEIN 
C  SHOEMAKER 

AMSRD  ARL  WM  T 
B  BURNS 

AMSRD  ARL  WM  TA 
W  BRUCHEY 
M  BURKINS 
W  GILLICH 
B  GOOCH 
T  HAVEL 
E  HORWATH 
M  NORMANDIA 
J  RUNYEON 
M  ZOLTOSKI 
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1  LTD 

R  MARTIN 
MERL 

TAMWORTH  RD 
HERTFORD  SGI 3  7DG 
UK 

1  SMC  SCOTLAND 
P  WLAY 
DERA  ROSYTH 
ROSYTH  ROYAL  DOCKYARD 
DUNFERMLINE  FIFE  KY  1 1  2XR 
UK 

1  CIVIL  AVIATION 
ADMIN  STRATION 
T  GOTTESMAN 
PO  BOX  8 

BEN  GURION  INTRNL  AIRPORT 

LOD  70150 

ISRAEL 

1  AEROSPATIALE 
S ANDRE 

A  BTE  CCRTE  MD132 
316  ROUTE  DE  BAYONNE 
TOULOUSE  31060 
FRANCE 

1  DRA  FORT  HALSTEAD 
P  N  JONES 

SEVEN  OAKS  KENT  TN  147BP 
UK 

1  SWISS  FEDERAL  ARMAMENTS 
WKS 
WLANZ 

ALLMENDSTRASSE  86 
3602  THUN 
SWITZERLAND 

1  DYNAMEC  RESEARCH  LAB 
AKE  PERSSON 
BOX  201 

SE  151  23  SODERTALJE 
SWEDEN 


1  ISRAEL  INST  OF  TECHLGY 
S  BODNER 

FACULTY  OF  MECHANICAL 

ENGR 

HAIFA  3200 

ISRAEL 

1  DSTO 

WEAPONS  SYSTEMS  DIVISION 
N  BURMAN  RLLWS 
SALISBURY 

SOUTH  AUSTRALIA  5108 
AUSTRALIA 

1  DEF  RES  ESTABLISHMENT 
VALCARTIER 
A  DUPUIS 

2459  BLVD  PIE  XI  NORTH 
VALCARTIER  QUEBEC 
CANADA 

PO  BOX  8800  COURCELETTE 
GOA  IRO  QUEBEC 
CANADA 

1  ECOLE  POLYTECH 
J  MANSON 
DMX  LTC 

CH  1015  LAUSANNE 
SWITZERLAND 

1  TNO  DEFENSE  RESEARCH 
R IJSSELSTEIN 
ACCOUNT  DIRECTOR 
R&D  ARMEE 

PO  BOX  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

2  FOA  NATL  DEFENSE  RESEARCH 
ESTAB 

DIR  DEPT  OF  WEAPONS  & 
PROTECTION 
B  JANZON 
R  HOLMLIN 
S  172  90  STOCKHOLM 
SWEDEN 


21 


NO.  OF 

COPIES  ORGANIZATION 


2  DEFENSE  TECH  &  PROC 
AGENCY  GROUND 
I  CREWTHER 
GENERAL  HERZOG  HAUS 
3602  THUN 
SWITZERLAND 

1  MINISTRY  OF  DEFENCE 

RAFAEL 

ARMAMENT  DEVELOPMENT 
AUTH 

M  MAYSELESS 
PO  BOX  2250 
HAIFA  3 1021 
ISRAEL 

1  TNO  DEFENSE  RESEARCH 
I H  PASMAN 
POSTBUS  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

1  B  HIRSCH 

TACHKEMONY  ST  6 
NETAMUA  42611 
ISRAEL 

1  DEUTSCHE  AEROSPACE  AG 

DYNAMICS  SYSTEMS 
M  HELD 
PO  BOX  1340 

D  86523  SCHROBENHAUSEN 
GERMANY 
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